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We investigate the anti-adiabatic limit of an anti- 
ferromagnetic S = 1/2 Heisenberg chain coupled to Einstein 
phonons. The flow equation method is used to decouple the 
spin and the phonon part of the Hamiltonian. In the effective 
spin model long range spin-spin interactions are generated. 
We determine the phase transition from a gapless state to a 
gapped (dimerized) phase, which occurs at a non-zero value 
of the spin-phonon coupling. In the effective phonon sector a 
phonon hardening is observed. 



Dedicated to Professor E. Miiller-Hartmann on the occa- 
sion of his 60th birthday. 



I. INTRODUCTION 

The interest in the spin-Peierls (SP) transition [1] has 
been renewed by the discovery of the first inorganic SP 
substance CuGe03 [2,3]. The SP instability of quasi one- 
dimensional spin systems is the analogon of the Peierls 
instability of quasi one-dimensional metals towards lat- 
tice modulations with a wave vector 2k p ■ Due to the 
coupling of the lattice degrees of freedom to the d = 1 
magnetic degrees of freedom the system can lower its en- 
ergy by undergoing a phase transition into a dimerized 
state. The loss in elastic energy is overcompensated by 
the gain in the magnetic energy. 

Most theoretical studies of SP transitions rely on an 
adiabatic treatment of the phonons (e.g. Ref. [4]). This 
means that the energy of the phonons, that are associ- 
ated with the lattice distortion, should be low compared 
with the spin-spin exchange coupling. In addition to the 
occurrence of a pre-existing soft phonon, one observes a 
phonon softening in the organic SP compounds. Fur- 
thermore, the phonon energy lo should be small com- 
pared with the resulting gap A in the spin excitation 
spectrum [5]. 

It has been recently emphasized [6] , that CuGeC>3 does 
not fulfill these conditions. Even more, a phonon harden- 
ing by about 5% — 6% is observed [7]. Uhrig [6] has devel- 
oped a flow equation approach [8] which is not founded 



on adiabaticity. The comparison with the DMRG data of 
Bursill et al. [9] has shown that the flow equation method 
works well, especially in the limit of small J/ui (cf. Fig. 1 
in Ref. [9]). In this paper we extend this approach. De- 
tails will be given below. 



II. FLOW EQUATION APPROACH 

The method of our choice is the flow equation approach 
introduced by Wegner [8]. This method is similar to 
Frohlich's approach [10], which projects a Hamilton op- 
erator to an effective Hamiltonian in one step applying 
an unitary transformation. Wegner modified the Frohlich 
transformation as he uses an infinite number of infinites- 
imal unitary transformations. 

The model we study reads 

H(£)=H s +Hsb+Hb (la) 
H s = ^(JiWS,S l+1 + J a (*)S,S, +a ) (lb) 

i 

H B = £ [u(£)b% + ^§-(b\b\ + bA)] (1c) 

i 

H SB = J2(MZ)bl + Ai(£)b i ) . (Id) 

i 

Here Si stands for the S— 1/2 spin operator on site i and 
bi (b\) destroys (creates) a phonon on site i. Hs is the 
Hamiltonian of a frustrated spin chain or J\-Ji model. 
If the frustration parameter a = J%j ' J\ exceeds a critical 
value of a = a c = 0.241167(5), then the model under- 
goes a quantum phase transition from a gapless state to 
a gapped phase [11-13]. We start with a Hamilton oper- 
ator (la) consisting of pure spin (Hs) and phonon (7Yb) 
parts and a spin-phonon coupling term 7Ysb and use flow 
equations to decouple the spin-phonon system. 

The continuous transformation is parameterized by a 
flow parameter £ £ [0,oo]. So 7i(i — 0) stands for the 
original Hamiltonian as given. We will end up with the 
effective Hamiltonian Ti.(£ = oo) which is more "diag- 
onal" in the way that the direct spin-phonon coupling 
has been rotated away. The infinitesimal generator rj(£) 
defines the unitary transformation via 



dH(£) 
d£ 



[V(£),H(£)] 
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A choice of r](£) proposed by Wegner is 
77W = [H d W,HW] , 



(2) 



(3) 
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where Tid is the part of the Hamiltonian which is re- 
garded as diagonal. By choosing H<x{£) = Hs(£) + Hb{£) 
and the "off-diagonal" interaction part as H d(£) = 
7~Lsb(£), our generator reads 77(f) = [Hd,H d\ = [Hs + 
, T~Csb] ■ In Chapter IV we present an alternative choice 
of 77, while the canonical version (3) is used throughout 
the main part of this paper. 
As initial conditions we choose 

Ji(0) = J and J 2 (0) = 0. (4) 

So the untransformed pure spin part is the normal 
Heisenberg Hamiltonian. As we will see later on, for 
non-zero values of the flow parameter £ long range spin- 
spin interactions arise. As a consequence, the next- 
nearest neighbour coupling constant J 2 in the effective 
spin model is finite. The phonon part Hb = H^ + H^ 
includes a non-boson- number conserving 

termW^. This 

choice is due to the fact, that the generator 77 = [Hs + 
H^\Hsb] would produce terms proportional to 
(b\b\) anyway. So for our purpose it is convenient to cm- 
bed this terms in H and renormalizc fi(£). Henceforth, 
we choose 

w(0) = cj = 1 and fi(0) = (5) 

and give all coupling constants in units of ujq. 

Hsb describes the coupling between spin and phonon 
system. The spin-phonon coupling mainly influences 
neighbouring sites. Hence the coupling operator Ai(0) 
typically consists of nearest neighbour spin products. An 
appropriate choice for Ai(0) is the difference coupling 

A(0) = g(SiS i+1 - SiSi_i) , (6) 

which does not contribute to the undistorted phase due 
to the translational invariance [6], i.e. (^(0)) vanishes. 
According to Ref. [8], A% should be normal-ordered, i.e. 
A(0) — > Ai(0) - (A(0)). This is guaranteed by the 
specific choice (6). By choosing different coupling oper- 
ators, various mechanisms of how the lattice distortions 
influence the exchange integral J can be investigated. 
In the case of the difference coupling, the exchange cou- 
pling directly depends on the "position" of a spin between 
the neighbouring ones. For CuGe03 this corresponds to 
changes of two neighbouring Cu-O-Cu binding angles, 
i. e. one of them is enlarged on expense of the other one. 

Sometimes a local coupling Wgg = g J2i SiSi+i(b\+bi) 
is used in connexion with CuGeO-3 (see e. g. Ref. [14- 
16]): Single harmonic degrees of freedom directly modify 
the magnetic interaction. In Ref. [16] the gap due to 
dimerization was investigated by means of a high order 
perturbation expansion. A flow equation approach to 
this coupling type will be presented elsewhere [17]. 

Applying the unitary transformation, not only the cou- 
pling constant in (6) changes (g — > a(£)). Additionally 
longer range and multiple site interactions arise. Their In- 
dependent coupling constants are denoted by b n {£). The 



aim of this approach is to disentangle phonons and spins. 
Hence in the limit t —* 00 the coupling operator Ai{£) 
vanishes, i. e. a(oo) = and 6„(oo) = 0. 

We introduce the Liouville operator £ for the commu- 
tation with Hs- £0 := [H$,0]. In extension to Ref. [6] 
the Liouville operator £,(£), the phonon frequencies ui(£) 
and n{£) are f dependent. We choose for the generator 77 

^)=^{[(£ + a;)A-Ml]6l + 

i 

[(£-w)^J + Mi]6i} (7) 

according to r\ = [Hs + Hb,Hsb}- The general flow equa- 
tion (2) leads to a flow equation for the coupling operator 
Ai (£) by investigating the terms linear in phonon opera- 
tors of [77(f), H{£)] 

^ =-{[(£ + c,) 2 - M 2 ] A- 2^4} ■ (8) 

Due to the complicated structure of Eq. (8) we restrict 
ourselves to a subspace of operators, which are allowed 
to appear in Ai{£). A systematic choice is to consider all 
terms arising in Ai(0) and £Ai(0). The following table 
summarizes the resulting set of operators. 



coefficient 


spin operator 


a{£) 




bo(£) 


iSi_i(Sj x Sj+i) 


bi(£) 


i[Sj(Sj+i x S l+2 ) + Sj_2(Sj_i x Si)] 


b 2 (£) 


i[Si(S i+ i x S i+3 ) + Si_ 3 (Si_i x Si)} 


bs(£) 


i{S t+2 (S t x S J _ 1 ) + S,_ 2 (S 4 x S i+ i)] 



In next order in £ the operators arising in £, 2 Ai(0) are 
included leading to 14 additional three spin operators 
of the type (Si x Sj)(Sfc x S;) and five additional long 
range two spin operators. Eq. (8) can now be divided in a 
set of differential equations for the spin-phonon coupling 
constants a{£) and b n (£). 

The truncation of the coupling operator A% in order n 
(by constructing the operator space with all terms ap- 
pearing in Ai(0), £Ai(0), £ 2 A(0), £ n Ai{0)) ne- 
glects effects of order n + 1 in J. We note here, that 
in this step no expansion in g occurs. As the contribu- 
tion to the flow of the neglected operators corresponding 
to higher order contributions in J is small compared to 
the one of the leading order terms, we expect the flow 
equations to converge for £ — > 00. This actually is the 
case except for large values of J/u>. 

For the phonon frequencies we find 

^ = (E([( £ + ^ A f ] + [A, (-£ + u)Al\) ) (9a) 

i 

^ = 2(£[£A-MUi])- (9b) 

i 

Only nearest and next-nearest neighbour spin products 
are kept and their mean field value, denoted by the outer 
brackets (•), is inserted. 
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The flow equations for J\ (£) and J2 (£) can be extracted 
from [r/,7Y SB ] 

E,{[(£ + ^A,4] +h.c.}(6+6) 
+ E i {[£A-Ml,A] +h.c.}(6+&t) (10) 
-E^A^^ + ^A-Mll+h.c.} . 

The quadratic boson terms are replaced by their expec- 
tation values. This is systematic in the sense of an ex- 
pansion in g. This mean- field step neglects fluctuation 
effects of order g 2 . As these terms only arise in g 2 -tcrms, 
the total error is of order g 4 (cf. Ref. [6]). We note here, 
that these expectation values are ^-dependent due to the 
^-dependence of u and fi. Finally, ^ is equal to 

the (next) nearest neighbour part of (10), respectively. 



1.2 




i 1 1 


1 _ 




1.0 


< — " " 




»(/) 




0.8 










0.6 










0.4 




J,(0)=0.6 
a(0)=g=0.3 






0.2 










0.0 


\ — 











0.0 5.0 10.0 15.0 20.0 

flow parameter / 



FIG. 1.: ^-dependent flow of the the coupling constants, the 
phonon frequencies and the coefficients in the ansatz for the cou- 
pling operator A in units of uiq ■ 



III. RESULTS 

Having collected all flow equations, the set of differen- 
tial equations (9 in order £ and 28 in order £ 2 ) have to 
be solved numerically. This is done via a standard fifth- 
order Runge-Kutta routine with adaptive step size con- 
trol. The change of the ^-dependent coupling constants 
for two succeeding Runge-Kutta steps is computed for 
all couplings. Convergence is assumed if the sum of the 
changes squared is less than 10~ 9 in units of ujq. Fig. 1 
shows for J/ujq = 0.6 and g/uJo = 0.3 the flow of all 9 
parameters. For intermediate values of I the additional 
four spin coupling terms arise (b n (£) 7^ 0). As £ be- 
comes larger, all spin-phonon coupling constants vanish, 
i. e. Hsb — » 0. Ji and J 2 reach their effective values 
in the limit I — > 00 and thus give rise to a frustration 
Qiefr = Ji S / m the decoupled spin model. This proce- 
dure is very intuitive in the sense that a flow diagram like 
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1 1 1 1 1 1 1 1 1 1— 
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FIG. 2.: Effective phonon frequencies w e ff (thin lines in the upper 
part), /i e fj (lower part) and A p h (thicker lines) for fixed values of 
the initial next-nearest neighbour coupling J = J\{i = 0). A p h = 

y/ oj^jj — /Ltj^j is the phonon frequency in the diagonalized version of 
the pure phonon part Hb- 

Fig. 1 shows, how the dressing of the spins with phonons 
induces a frustration a C R > 0. 

As a by-product, we find an increased phonon fre- 
quency in the effective phonon part of the model. This 
phonon hardening can be interpreted qualitatively as a 
level repulsion between the high energy phonon system 
and the low energy spin system. Recently, it has been 
pointed out by Gros and Werner [18] that this feature is 
consistent with the RPA-approach by Cross and Fisher 
to the spin-Peierls transition. Fig. 2 shows for four values 
of J the effective phonon frequencies in dependence of g. 
As an appropriate measure of the phonon hardening, the 
frequency A pri = yj w 2 ff — of the diagonalized version 

— > J2i ^phfljpi + A ph ,o is depicted. 

The phase diagram of the spin-Peierls anti-ferromag- 
netic chain of spins interacting with phonons has been 
investigated by Bursill et al. via DMRG [9], by Uhrig 
via flow equations in leading order in g/oj and J/u> [19], 
by Weifie et al. via fourth order perturbation theory and 
Lang-Firsov transformation [20] and by Sun et al. via 
mean field and renormalization group methods [21]. In 
Fig. 3 results of various flow equation approaches are 
shown. In Fig. 4, we compare our results to those of the 
other approaches. 

Our calculation is based on the idea that the transition 
into the ordered phase does not occur due to the soften- 
ing of a phonon but due to the tendency of the effective 
spin model towards dimerization. Phonon-induced frus- 
tration above its critical value a e s > a c drives the dimer- 
ization. The phase transition line is therefore determined 
by solving 

a c - ( lim a) I =0 (11) 
for fixed values of J with respect to g. In this way we 
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FIG. 3.: Zero temperature phase diagram of the spin-Peierls anti- 
fcrromagnctic chain of spins interacting with phonons. For small 
values of the spin-phonon coupling g/u>o the system is gapless. For 
large g/u>o the system is dimerized and has an energy gap. The 
curves denoted with ct sw and d sw are the results of two modified 
approaches and will be discussed in Chapter IV. 

find the critical spin-phonon coupling g c in dependence 
of the nearest neighbour spin coupling J. For values of 
J/u up to 0.35 the extended flow equation approach im- 
proves the results from Ref. [6]. For intermediate values 
of J/u the curve crosses the DMRG result and predicts 
critical couplings, which must be considered as too low. 
It is interesting that the result of Ref. [6] gives a good 
description of g c up to J/u w 10 although it is system- 
atic only to order g 2 /u 2 and g 2 J/u 3 . We presume that 
the two leading orders comprise already the main compe- 
tition between dimerization and delocalization. Further 
orders start to build in resonance effects. But they do not 
succeed for all J/u because the effective interactions are 
restricted to small distances only in our ansatz. It is not 
uncommon (cf. mean field treatments) that the leading 
orders yield physically reasonable though not accurate 
results whereas systematic improvements are difficult to 
implement (see Ref. [22] for an expansion in the inverse 
coordination number as example). 

The inclusion of the next order in the systematic ex- 
pansion by WeiBe et al. leads to slight improvements for 
J/u < 0.15 which is an estimate of the convergence ra- 
dius of the series expansion. But the higher order terms 
produce too large critical couplings even for intermediate 
values of J/u [20]. Furthermore, WeiBe et al. used the 
Lang-Firsov transformation to combine adiabatic and an- 
tiadiabatic effects in their theory. This leads actually to 
an additional increase of the critical couplings for inter- 
mediate J/u (cf. Fig. 6 in Ref. [15]). We conclude that 
the extended flow equation approach used here has im- 
portant advantages over the systematic series expansion. 
This is due to its renormalizing property, i. e. the feed- 
back of the ^-dependence of all couplings. In this way 
the region of validity is increased up to J/u « 0.35. 



FIG. 4.: Phase dia gram like Fig. 3. Various previous results are 
compared with the phase transition lines obtained via flow equa- 
tions or DMRG. 

For values of J/u < 1 Sun et al. use g c /u = \J J/(2u) 
[21,23]. This is smaller by about a factor l/\/2 than the 
first order result obtained by Kuboki and Fukuyama [24] , 
who find g c /u — y / a c /(l/2 — a c ) ■ J/u w y/ J/u. As 
can be seen from Figure 4, the results of Sun et al. do 
not compare so well with those of the other approaches. 
For J/u > 2 they find critical couplings g c /u s=s 0.6, 
which are about twice as large as the g c -values computed 
by DMRG [9]. For increasing J/u, Sun et al. predict 
decreasing g c above J/u w 5 (not shown) in contrast to 
the DMRG data. 

While for not too large values of J/u the extension of 
the flow equation approach to higher orders in g/u pro- 
vides a fast method to retrieve data as precise as DMRG 
calculations, the large J/u regime cannot be described 
so easily. At J/u ~ 0.77 the effective Hamiltonian is 
no longer hermitian and the description breaks down. 
The next systematic extension is to enlarge the operator 
subspacc for Ai(l) to order £ 2 . The result once again 
improves the calculations in order £ for small J/u, but 
cannot circumvent the "downturn" at J/u w 0.4. We 
will argue below that the inclusion of longer range inter- 
actions is necessary for J/u > 0.35. 

IV. MODIFICATIONS 

To enlarge the region of validity of the flow equation 
approach to larger Js, we introduce some modifications. 
In a first attempt the restriction to a given operator sub- 
space can be dropped for intermediate terms arising in 
(8-10) via ZAi- This circumvents the downturn of the 
phase transition line but leads to too large critical cou- 
plings g c (J). As these modifications changed the solu- 
tions of the differential equations drastically the influence 
of longer range interactions present in these intermediate 
terms can be judged as important. To include them in 
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a systematic way two additional couplings J3 and J4 are 
introduced. They correspond to spin interactions SiS^+3 
and SiS i+ 4. The flow equations for Js(£) and J±{£) can 
be extracted from (10). As the critical coupling a c is only 
known for the Ji-J 2 -model we combine these four cou- 
plings to two effective couplings Jf(£) and J? w 



thus ■ 



'(£) 



(£) and Jl w (£) and 
J| w / Jf w . This is done by using linear spin 



wave theory [25]. The effective couplings Jf w and J| w 
are chosen such that the spin wave dispersion with Jf w 
and J| w equals the one for Ji, J 2 , J3 and J 4 . Thus we 
obtain 



jr 



TSW 

J 2 — 



TSW 

J 2 



9J 3 )(Ji 

-4J 4 



Ji 

W 



9J 3 

J 2 - 



f J3) 
4J 4 ) 



9J 3 



(12a) 
(12b) 

(12c) 



The curve of Fig. 3 denoted with a sw is the solution of 
= 0. The downturn does not occur any more, 



a c -a"„ 



but the critical couplings are too large. As and =4r 
are solved separately from all other coupling constants, 
i. e. J3 and J 4 do not couple to the other differential 
equations, another method denoted with d sw is intro- 
duccd. Differentiating (12a) and (12b) one finds —A- 



and 

dj„ 



d£ 



These two derivatives do depend on all J„ and 
-gj- (n — 1, . . . ,4). They are calculated in each Rungc- 
Kutta step. This provides a coupling of J3 and J4 to the 
other flow equations. The results from this procedure 
show again the downturn at J/luq s=s 0.4 and produce 
worse data for small J/ uu than the unmodified approach. 

The inclusion of long range interaction terms modifies 
the results drastically. So they are important and have 
to be treated correctly. It seems that one has to extend 
the flow equation approach to very high orders in the 
Liouvillc operator £, i.e. the subspace for the coupling 
operator Ai{£) should be fairly large, if one wants ac- 
curate data for large values of J/oj. For J/uj < 0.35, 
however, the method presented works very accurately. 

Due to the complicated spin commutators and prod- 
ucts in the flow equations an extension to higher or- 
ders in £ is cumbersome. A further hint of how this 
could be achieved more easily shall be given. Recently 
non-standard choices of the generator rj have been used 
successfully [26,27]. The main point is to introduce a 
quantum number Q and define the generator via r](£) — 
[Q, H(£)}. Thus [Hes, Q] = 0, i. e. Q is conserved with re- 
spect to H c ff. For our problem we introduce the phonon 
number Q = J^bJ&i. With this choice the generator 
reads r\ = ^2AAib\ — A\bi) and does not contain the 
Liouville operator. The flow equations (for T = 0) are 
pretty simple 



^i = -(£ + c)A 
du> n ,. , 

^ = 2<[A,A f ]> 

[r 7 ,W S B]--2^A t A 



(13a) 

(13b) 
(13c) 



We investigated these flow equations in order £ and £ 2 . 
The results are similar to those of the standard choice of 
rj. Although the standard choice (3) leads to slightly bet- 
ter results for small couplings J/uj, the alternative gener- 
ator could be used to investigate very large subspaces for 
Ai(£), as only £ enters linearly in (13a) and it is absent 
in (13b,13c). 



V. SUMMARY 

We used the flow equation method to rotate the spin- 
phonon coupling of a spin-Peierls Hamiltonian away. In 
this way we derived effective spin and phonon models. 
The spin-phonon coupling induces longer range inter- 
actions, i. e. produces a non-vanishing frustration, and 
renormalizes the original couplings. The effective phonon 
frequency is enlarged. This corresponds to the phonon 
hardening observed for CuGc03. The phase diagram de- 
scribing the transition from a gapped to an ungapped 
state is compared with DMRG results. For large values 
of J/uo the results deviate. In this regime very long range 
interactions appear to be relevant and the transformation 
of the spin-Peierls chain onto a rather local effective spin 
model must be considered as not sufficient. The flow 
equation and the DMRG results coincide for J/u> < 0.35. 
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